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JL.he long-standing problem of equatorial superrotation in a slowly rotating planetary 
atmosphere has been solved recently. This article presents a systematic, comprehensive, 
and analytical examination of the momentum budget that maintains stable superrota- 
tional winds in a given planetary atmosphere. After reviewing general approaches used in 
modeling comparative planetary atmospheres, I describe the major dynamical processes 
that maintain stable equatorial superrotation for the atmospheres of Venus, Saturn’s 
moon Titan, and Earth. Under appropriate conditions, all three bodies could have equa¬ 
torial superrotational winds greater than 100 ms -1 . Venus’ equatorial superrotation of 
118ms” 1 is maintained at its cloud-top level at =65 km. My model predicts that, to main¬ 
tain an equatorial superrotational jet stream in Titan’s stratosphere, its main haze layer, 
which absorbs solar radiation, would need to be centered at an altitude of 288 km (corre¬ 
sponding to 10 Pa pressure). It is hoped that ongoing observations by the Cassini-Huygens 
mission to Saturn will be able to validate this model prediction. Strong equatorial superro¬ 
tational winds could also be induced on Earth by optically thick dust clouds, such as those 
arising from an asteroid impact or extreme volcanic event. Superrotation would greatly 
enhance meridional transport in the Earth’s atmosphere, which could have accelerated 
the globalization of an environment unfavorable to the survival of the dinosaurs. 


THE SCIENCE OF COMPARATIVE 
PLANETARY ATMOSPHERES 

This article gives a brief overview of modeling stud¬ 
ies of planetary atmospheric physics, focusing mainly 
on atmospheric motions associated with dynamics, 
and reports significant progress in solving the long¬ 
standing dynamical problem of equatorial superrotation 


in a slowly rotating planetary atmosphere. 1 The discus¬ 
sions are presented from the perspective of comparative 
planetary atmospheres rather than for a specific planet. 
The scientific merits of conducting studies of compara¬ 
tive planetary atmospheres are illustrated based on both 
general rationales and specific examples. 
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Over the past few decades, our understanding of the 
planets (including planetary bodies such as Saturn’s 
moon Titan) and their atmospheres has advanced from 
a highly subjective art to a purely scientific process, sup¬ 
ported by a large number of measurements made by 
various space missions. Our current knowledge indi¬ 
cates that while planetary atmospheres share physical 
processes, they exhibit diverse characteristics. The goal 
of modeling comparative planetary atmospheres is to 
understand a variety of atmospheric phenomena gov¬ 
erned by common physical laws. 

One way to illustrate the importance of compara¬ 
tive planetary atmospheres is to consider the Earth’s 
atmosphere from first principles. There are two clear- 
cut methodologies in the physical sciences: (1) theo¬ 
retical modeling based on observations and (2) labo¬ 
ratory experiments. The first method starts from a set 
of fundamental physical laws that govern the states 
and motions of the atmosphere. These physical laws 
are expressed mathematically and are often analyzed 
and solved numerically. Solutions are compared with 
measurements, from which our understanding of the 
Earth’s atmosphere can be steadily improved. The 
essence of the second methodology is the simulation 
of processes in the laboratory, where one can con¬ 
trol and change the external conditions that govern 
the behavior of the particular physical system under 
study. One must therefore be able to reproduce an 
approximate physical model of the real system in the 
laboratory where it can be controlled—a require¬ 
ment that clearly cannot be met for the Earth’s entire 
atmosphere. In fact, atmospheric research is mostly 
an observational science rather than an experimental 
science. 

With an observational science, the critical process 
of hypothesis testing is severely curtailed because of 
our inability to control the external variables in labora¬ 
tory experiments designed to simulate atmospheric pro¬ 
cesses. However, planetary atmospheres with remark¬ 
ably different astronomical conditions, which result 
in significantly different atmospheres in terms of their 
physical and chemical properties, serve as ideal “labo¬ 
ratories.” They provide a natural set of experiments, 
allowing us to test our basic understanding of physical 
processes that are governed by the same set of physical 
laws, yet under different external conditions. Therefore, 
comparing and contrasting different atmospheres in the 
solar system enables us to test our basic understanding 
of atmospheric processes in general and can also help us 
gain insight into the Earth’s atmosphere. For example, 
the well-known phenomenon of the Antarctic “ozone 
hole” was found to be caused by the catalytic destruc¬ 
tion of ozone by chlorine and other halogen-contain¬ 
ing components. This discovery benefited greatly from 
theoretical studies of the photochemistry of HC1 in the 
atmosphere of Venus. 2 


The Earth’s atmosphere has been studied extensively, 
mainly by meteorologists. Yet meteorologists have dif¬ 
ferent objectives and approaches from those studying 
comparative planetary atmospheres. The main objective 
of studying Earth’s atmosphere is to make better predic¬ 
tions of its future state. In weather forecasting as well as 
climate prediction, one likes to know the physical state 
of the atmosphere at any given time. The prediction 
problem is often treated as an “initial value problem” 
in mathematical terms. So-called data-model assimila¬ 
tion, which appropriately incorporates current measure¬ 
ments into a continuously running numerical model, 
often plays an important role in prediction. Conversely, 
in comparative planetary atmospheres, one often cares 
more about the physical mechanisms that explain or pre¬ 
dict atmospheric phenomena. In these cases of “mecha¬ 
nistic studies,” the problem is often treated as a so-called 
forcing-response problem, and data-model assimilation 
usually does not enter directly into the formulation and 
solution procedures in numerical modeling. The physical 
mechanisms in comparative planetary atmospheres usu¬ 
ally include those that explain the causes of atmospheric 
composition, origin, evolution, and dynamics. 

This article focuses mainly on modeling studies of 
atmospheric motions associated with dynamics and 
reports on recent progress in solving the long-standing 
dynamical problem of the equatorial atmosphere rotat¬ 
ing much faster than its solid planetary body (superrota¬ 
tion) in a slowly rotating planetary atmosphere. 

OVERVIEW OF ATMOSPHERIC 
MODELING STUDIES 

In modeling studies of atmospheric states and 
motions, significant advances can be made by investigat¬ 
ing common atmospheric processes, such as energy bal¬ 
ance, the momentum budget, the interaction between 
waves and the mean flow, etc., as they apply to different 
planetary settings. Atmospheric features observed on 
different planets can frequently be explained by under¬ 
standing common processes governed by common physi¬ 
cal laws. 

The basic physical laws describing the conservation of 
momentum, energy, and mass are the same for all plane¬ 
tary atmospheres. The mathematical expressions of these 
laws are called “primitive equations” (see Box 1, “Invari¬ 
ant Primitive Equations”). Often, one does not solve the 
complete set of primitive equations directly when study¬ 
ing atmospheric motions. Box 1 gives one example of 
deriving an important feature of planetary atmospheric 
motions by a set of simple sign transformations. Numeri¬ 
cal models for other planetary atmospheres are usually 
derived from those for the Earth’s atmosphere. An invari¬ 
ant characteristic demonstrated in Box 1 indicates that it 
is the magnitude rather than the direction of planetary 
rotation that matters. 
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BOX 1: INVARIANT PRIMITIVE EQUATIONS UNDER A COORDINATE TRANSFORMATION 

The so-called atmospheric primitive equations for the V = velocity vector (u, v, w), 

Earth’s atmosphere are expressed in spherical coordinates that x, y = eastward and northward distance, respectively, 
describe the fundamental conservation laws for atmospheric Z= altitude, 

momentum, energy, and mass, 3 that is, p = air density, 

<t>= latitude, and 

momentum equations Q = Earth > s rotation rate 

The terms F rx , F ry , and F r? in Eqs. 1—3 and Q in Eq. 4 rep¬ 
resent the net forcing to the momentum and internal energy, 
respectively. The total derivative D/Dt with respect to time 
t and wind divergence V*V in spherical coordinates that 
adopted appropriate approximations for the Earth’s atmo¬ 
sphere are defined as follows: 
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the thermodynamic energy equation 


D(c T) p ^ 

-U«-A-c£vv = Q; 


Dt 


(4) 


and the continuity equation 
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In these equations, 

a= Earth radius, 

c v = specific heat at constant volume, 
c v T = internal energy, 
g= gravity, 
p= air pressure, 

T = air temperature, 
t= time, 

u, v, <w= eastward, northward, and upward velocity, 
respectively, 


„ TT du 1 d(cos </>v) die ... , 

V-V =-1- 1 —H-. (6b) 

dx cos (f> dy dz 

Equations 1—6 can also be applied to studying the atmo¬ 
spheric motions of other planetary bodies. The inclination of 
the Earth’s equator to its ecliptic is 23.45° and its rotational 
direction is from west to east. Some planets or planetary bodies 
have negative inclinations. For example, Venus’ inclination 
is -2.7° and Pluto’s is -57.54°. Venus and Pluto rotate retro¬ 
grade with respect to Earth: from east to west. 

It can be easily verified that Eqs. 1-5 are invariant after 
making the following sign transforms: Q. —h - Q, x —» — x, and 
u —h — u; i.e., if our Earth were rotating from east to west, the 
magnitude of its atmospheric motion would remain the same 
except that its wind direction would be reversed. In other 
words, the atmospheric states as solutions of a rotating dynam¬ 
ical system can be more effectively described by a coordinate 
system based on a rotating axis rather than fixed east-west 
coordinates. Therefore, when studying atmospheric motions 
from the standpoint of comparative planetary atmospheres, 
one can define prograde (retrograde) wind to be motions in 
the same (opposite) direction as the rotational direction of 
a planet. Likewise, prograde (retrograde) acceleration corre¬ 
sponds to a momentum source (sink) that generates (destroys) 
the east-to-west winds in the atmosphere of Venus and the 
west-to-east winds in the atmospheres of Earth and Titan. 


In terms of the number of space variables explicitly 
included, the state, structure, and circulation of a plan¬ 
etary atmosphere can be studied by 0-, 1-, 2-, or 3-D 
models. As the number of model dimensions increases, 
more realistic physical processes that are often parame¬ 
terized (i.e., the detailed processes are expressed in terms 
of simplified empirical formulas with a few hulk param¬ 
eters) in lower-dimensional models are included. Even 
though more physical processes are explicitly modeled 
in a higher-dimensional model, the number of adjust¬ 
able and empirical model parameters still increases for a 
fixed number of observations because of the much faster 
increase in the number of model gridpoints. There¬ 
fore, the computational resource requirements increase 


dramatically as one proceeds to higher-dimensional 
models. 

There is always a fundamental trade-off between 
model efficiency and model robustness in terms of 
choosing the appropriate model dimensionality and the 
appropriate physical processes to he explicitly included 
in the model. Three-dimensional models are the most 
complicated, addressing all interactions among dynam¬ 
ics, photochemistry, radiation, and transport. However, 
the merit and success of using 3-D models for studying 
comparative planetary atmospheres depend largely on 
the availability of data that can be used to initialize and 
constrain the model; that is, the availability of data pro¬ 
vides constraints on the complexity of models. Here, I 
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give a brief description of atmospheric models in lower 
dimensions. 

Zero-dimensional models describe atmospheric states 
that are independent of spatial variables (i.e., location). 
A simple example of a 0-D model is the calculation of 
the effective temperature T e of a planet, which is defined 
by an energy balance between the incoming solar flux 
and outgoing thermal radiation 4 : 


TTtl (1 — Cl )S 


-= 4'ITCl 2 0'T 4 


( 1 ) 


where 

a = planet radius (6371 km for Earth), 
a p = planetary albedo (0.30 for Earth), 
d = Sun-planet distance (1.00 AU for Earth), 

S = solar constant at 1 AU (1366 W m _2 )> 

T g = effective temperature, and 

(X = Stefan-Boltzmann constant (5.67 X 10 _ ® W rrD 2 ). 

The left-hand side of Eq. 1 denotes the solar energy 
intercepted and absorbed by the planet disk, whereas 
the right-hand side denotes the thermal energy radiated 
back out to space by the planet as a blackbody radiator. 
Equation 1 can be considered a degenerate form of Eq. 
4 in Box 1: Q = 0. 

Table 1 lists several values of the effective tempera¬ 
ture T e and the mean surface temperature T s for a few 
selected planets. Because of the greenhouse effect, T s 
is often greater than T. “Greenhouse effect” refers to 
the planet’s surface receiving both direct shortwave 
solar radiation that transparently passes through the air, 
which would yield T , and longwave IR radiation emit¬ 
ted from the air that absorbs and thus prevents the IR 
radiation from below from returning to space, leading 
to T s > T. The higher T s than T e for Jupiter is mainly 
due to its internal heat source. In addition to Jupiter, 
Saturn and Neptune also possess substantial internal 
heat sources compared to the solar energy they receive. 
For a planet with a very thin atmosphere such as Pluto, 
T e is approximately the same as the planetary T. For 
Pluto’s icy surface, I assume a relatively large value of 
albedo: a = 0.67. On the basis of the above energy bal¬ 
ance equation, it can be found that T e = 38.7 K, which 
is consistent with the value derived from coincident 
measurements . 5 

One-dimensional models often refer to radiative-con¬ 
vective or radiative-conductive models that determine 
the vertical (z) thermal structure or photochemical 
(box) models that determine the vertical distribution 
of chemicals in the atmosphere. Atmospheric radiative 
transfer is key to understanding and developing various 
1-D models of the vertical thermal structure of plan¬ 
etary atmospheres . 6 Measured vertical profiles of tem¬ 
perature, air density, and chemical composition provide 


Tablet. Effective temperatures T e and mean surface 
temperatures T s for select planets. 

Sun-planet 



distance, 

VAU) 

Albedo, 

a P 

T e (K) 

T,(K) 

Venus 

0.72 

0.77 

227.0 

750.0 

Earth 

1.00 

0.30 

255.0 

280.0 

Mars 

1.52 

0.15 

217.0 

225.0 

Jupiter 

5.20 

0.58 

98.0 

138.0 3 

Titan 

9.57 

0.30 

82.0 

94.0 

Pluto 

29.70 b 

0.67 

38.7 

38.7 


a At cloud top (pressure = 1 bar or 10 5 Pa). 
b In May 1993. 


constraints on atmospheric energy exchange, photo¬ 
chemistry, and transport within the atmosphere and 
possibly also on thermodynamic and chemical states at 
the lower boundary or in the deep atmosphere. Since 
the distribution of minor photochemical species often 
exhibits the largest gradient in altitude, 1-D models are 
useful tools for interpreting the measurements of minor 
chemical constituents. Note that the major components 
of planetary atmospheres are generally uniformly dis¬ 
tributed globally, such as N 2 for Earth and Titan, C0 2 
for Venus and Mars, and H 2 for giant planets. Therefore, 
their distribution can often be derived by 0-D models 
alone on the basis of thermodynamics and from the 
perspective of solar system evolution . 2 Since the minor 
constituents of planetary atmospheres (e.g., O 3 for Earth 
and Mars, HC1 for Venus, and CH^ for giant planets) 
vary the greatest in the vertical direction and are sensi¬ 
tive to the vertical thermal and dynamical structure of 
the atmosphere, a 1-D model is the minimum require¬ 
ment for determining the vertical profiles of minor spe¬ 
cies in photochemical-diffusive equilibrium in a plan¬ 
etary atmosphere. Various 1-D models for different plan¬ 
ets are presented in a recent monograph by Yung and 
DeMore . 7 

Two-dimensional models, which contain two spatial 
variables (altitude and latitude), have been useful tools 
for mechanistic studies of the atmospheres of Earth and 
other planets. For a rotating planetary atmosphere, zonal 
mean (latitude bands averaged across all longitudes) 2-D 
models can be considered a good trade-off, mainly 
because large-scale atmospheric structure and circulation 
vary the least along latitude bands. In a zonal mean 2-D 
model, atmospheric state variables such as temperature, 
winds, and chemical species distributions are expressed 
as functions of time t, latitude <J>, and altitude z- Cur¬ 
rent computational capacity allows one to couple com¬ 
plicated photochemistry with dynamics through trans¬ 
port in a 2-D model. However, this kind of model can 
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only be well tested for the Earth’s atmosphere because of 
the availability of databases of dynamical state variables 
and chemical species. The most observable and stable 
dynamical features for planetary atmospheres are zonal 
jets at various altitudes and latitudes. Important ques¬ 
tions in planetary atmospheric dynamics remain: How 
is an atmospheric jet stream maintained on a particular 
planet? Why is the jet stream for that planet located at 
a particular latitude and altitude? What determines the 
magnitude of the maximum wind at the jet center? Is 
the observed jet stream a stable, steady physical state or 
a transient realization of an oscillatory state? 

MAINTENANCE OF EQUATORIAL 
SUPERROTATION IN A PLANETARY 
ATMOSPHERE 

A long-standing problem in planetary atmospheric 
dynamics is the maintenance of equatorial super¬ 
rotation in the atmospheres of slowly rotating plan¬ 
ets. Strong equatorial superrotational zonal winds of 
=100 m s _1 (Fig. 1) were first observed at Venus’ cloud- 
top level in the 1960s, 8 consistent with an atmospheric 
angular velocity nearly 60 times greater than the rota¬ 
tion of the solid planet! Several mechanisms have been 
proposed to explain the generation or maintenance of 
the superrotation of clouds on Venus. 



Figure 1 . A UV image of thick clouds on Venus taken by the Pio¬ 
neer Venus Orbiter (26 February 1979). Venus rotates very slowly, 
taking 243 days to complete one revolution. By tracking the cloud 
features, it is found that the clouds travel quickly around the planet, 
taking only 4 days to circulate once, indicating superrotational 
winds at the cloud level. Both the winds and the solid planet rotate 
from east to west, opposite to the direction on Earth. The Y-shaped 
cloud patterns move from right to left in the figure. 


I have already shown in Box 1 the primitive equa¬ 
tions for atmospheric motion applicable to all planets. 
To examine the dynamics of a zonal jet stream, one 
needs to take a zonal average over Eqs. 1-5 in Box 1 
and evaluate the zonal momentum budget of all terms. 
Since the zonal momentum equation describing the 
zonal mean winds contains 10 terms (see Eq. 2 below), a 
rigorous and quantitative study of the superrotation has 
always been difficult. Most theories have qualitatively 
or quantitatively described possible mechanisms arising 
from one or two particular dynamical processes corre¬ 
sponding to one or two terms in the zonal momentum 
equation. Also, most theories have studied equatorial 
superrotation based on a set of atmospheric parameters 
specific to Venus. Two well-known theories that explain 
the maintenance of Venus’ equatorial superrotation are 
(1) meridional momentum transport from mid-latitudes 
by eddy mixing 9,10 and (2) momentum pumping by ther¬ 
mal tides. 11,12 Careful examination of the existing theo¬ 
ries shows that they either have difficulties leading to 
a contradiction with fundamental physical laws or are 
incomplete, leading to incorrect predictions of the mag¬ 
nitude of the superrotating winds. 

To rigorously solve problems of large-scale circula¬ 
tion, such as equatorial superrotation in a planetary 
atmosphere, a 2-D model is necessary to describe the 
zonally averaged winds. By taking a zonal average of the 
flux form of Eq. 1 in Box 1,1 arrive at the following mean 
zonal momentum equation 1 : 


I 2 

pu t =--- ipv'u' COS ~ (pu/u')^ — pvou~ 

acos A <p 

(I) (II) (III) (IV) 

— 2£2 cos <f>{pw) — a R pu + 2X2sin^>(pv) 

(V) (VI) (VII) (2) 

(wCOS(/>)^ _ p W ' u ' puiu 

_ pv , 

a cos cp a a 

(VIII) (IX) (X) 

where the overbars and primes denote the zonal aver¬ 
ages and deviations from them, respectively. I have also 
replaced the momentum drag term F by a Rayleigh 
frictional drag term —a R pu with the Rayleigh frictional 
coefficient a R . Equation 2 shows that the mean zonal 
momentum equation contains 10 terms (I-X), each rep¬ 
resenting a different dynamical process contributing to 
the momentum sources and sinks. These terms have 
never been systematically and self-consistently evaluated 
in the same dynamical frame. Therefore, the exact phys¬ 
ical mechanism that maintains equatorial superrotation 
on Venus has remained a mystery, and the previously 
proposed hypotheses are mutually inconsistent. This 
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Figure 2. (a) Oscillating heating (+) and cooling (-) centers move opposite to the direc¬ 
tion of the zonal wind. The phase surfaces (tilted lines) and air parcel trajectories (arrows) 
of the excited tidal waves cause a negative mean momentum flux above the heating layer 
and a positive one below the heating layer. Vertical momentum flux convergence across 
the heating layer produces an acceleration of the mean zonal flow, (b) Momentum flux 
convergence within the heating layer can also be represented as the vorticity flux. A posi¬ 
tive correlation between w' and also leads to the acceleration of the mean zonal flow, 
(c) The incoming branches of the Hadley circulation cell cause the convergence of the 
momentum-depleted air mass from the outside into the region of the superrotational jet, 
whereas the outgoing branches cause the divergence of the momentum surplus air mass 
from the jet center to the region outside. The Hadley circulation is a zonally averaged 
thermal circulation that consists of warm air rising at the equator, moving poleward, then 
descending at the high latitudes and moving equatorward at the lower altitudes. 


lack of resolution of the superrota- (a) a 
tion issue on Venus has also caused 
further confusion and debate con¬ 
cerning equatorial superrotation in 
Titan’s stratosphere (e.g., Refs. 13- 
15). Titan is also a slowly rotating 
planetary body, which, similar to 
Venus’ cloud layer, has a haze layer 
that absorbs a significant portion of 
its incident solar radiation. 

Recently, I solved the problem of 
equatorial superrotation in a plan¬ 
etary atmosphere by specifying the 
2-D zonal winds in an analytical 
form that resembles a jet stream. 1 
As a result, each of the 10 terms in 
Eq. 2 can be analyzed systematically, 
first by scale analysis, then through 
analytical formulations and numer¬ 
ical evaluations, and finally by 
comparison with well-known plan¬ 
etary atmospheres. Key parameters 
of the analytical formulations are 
consistent with fundamental physi¬ 
cal concepts and agree with obser¬ 
vational constraints. Furthermore, 
this theory has been developed from 
the general perspective of equatorial 
superrotation for any given plane¬ 
tary atmosphere. The analytical 
forms of the approach explicitly 
show how an equatorial superrota¬ 
tion depends on various external 
and internal parameters. This general theory has been 
applied to the atmospheres of Venus, Titan, and Earth. 

For the atmospheres of these three bodies, which 
have or could have equatorial superrotation, the terms 
that make the major contributions to the momentum 
budget in Eq. 2 are III, VI, and VIII. Dynamically, 
HI = — (pw’u) represents the convergence of the verti¬ 
cal eddy momentum flux caused by tidal wave pump¬ 
ing, i.e., the dynamical effect of waves on mean flow 
through the nonlinear interaction between the two. 
The physical mechanism is shown schematically in 
Fig. 2a. When an oscillatory heating source caused by 
the absorption of solar radiation moves in a direction 
opposite to the zonal flow, the phase surfaces (tilted 
lines) and air parcel trajectories (arrows) of the excited 
tidal waves slant differently above and below the heat¬ 
ing layer. As a result, a vertical momentum flux conver¬ 
gence arises across the heating layer, which produces 
an acceleration of the mean zonal flow. Feovy 12 also 
noted that the momentum flux convergence can be 
rewritten as 


- (pw'u') = pw'g ', (3) 


where = — u' + «/ is the horizontal vorticity (along 
the y direction) of the perturbation field that corre¬ 
sponds to tides in this study. Equation 3 shows that the 
momentum flux convergence within the heating layer 
can also he represented as a vorticity flux. Figure 2b illus¬ 
trates how a positive correlation between w' and 
leads to an acceleration of the mean zonal flow at the 
center as the vortex parcels move in an organized way. 
Though the schematic depictions of III shown in Figs. 2a 
and 2b are apparently different, they are actually identi¬ 
cal in terms of mathematics and physics. 

Since the momentum source III produces prograde 
acceleration of the equatorial superrotational winds, 
terms VI and VIII produce retrograde acceleration to 
balance III. The momentum sink VI = — cr R pu due to 
parameterized Rayleigh friction represents the effect 
of small-scale eddies. These eddies cascade the kinetic 
energy of the flow from large-scale motions to smaller 
ones until completely dissipated by the molecular viscos¬ 
ity. The momentum sink VIII = — pv(acos(/>) -1 (ucosd>)^ 
represents the advection of horizontal wind shear by 
the horizontal branches of the Hadley circulation cell. 
The dynamical mechanism of momentum sink VIII 
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Figure 3. Dimensionless relationship between the wave pumping and momentum sinks 
as a function of fhe magnitude of the zonal wind at the superrotational jet center. The 
red curve represents the momentum source caused by wave pumping. Three straight 
lines represent relative momentum sinks under different parameter settings. A steady 
and stable superrotational solution is possible when two intersections exist between the 
source curve and a sink line. 


decelerates the wind at the jet 
center as shown schematically in 
Fig. 2c. The incoming branches of 
the Hadley circulation cause the 
convergence of the momentum- 
depleted air mass into the region of 
the superrotational jet, whereas the 
outgoing branches cause the diver¬ 
gence of the momentum surplus air 
mass from the jet center to the out¬ 
side region. 

As a result of wave/mean-flow 
interaction, the functional depen¬ 
dence of the wave pumping III on 
the zonal wind at the jet center u Q 
is highly nonlinear, as shown sche¬ 
matically in Fig. 3 by the red curve 
(see Box 2, “Stability Analysis”), 
which vanishes when u Q gets too 
small or too large. However, both 
VI and VIII are linearly propor¬ 
tional to u Q , so the sum of the two 
is a straight line passing through 
the origin. A steady and stable solu¬ 
tion exists where the momentum 
sink line intersects the momentum 
source curve (Fig. 3). The peak 
value of the red curve and slope of 
the straight lines in Fig. 3 depend 
on the strength of the solar heating 
rate Q, since both the amplitude 
of the tidal waves ( u\ v', w') and 
the strength of the Hadley circula¬ 
tion (v, w) are linearly proportional to Q. However, III 
depends quadratically on the wave amplitude, whereas 
VIII depends linearly on v. Hence, the relative strength 
of the momentum sink VIII is inversely proportional to 
Q. Note that VI is independent of Q. Thus, the slope of 
the straight lines in Fig. 3, which measures the relative 
strength of the total momentum sinks, increases with 
increasing d sp , since less solar energy is received when a 
planet gets farther away from the Sun. 

I show in Table 2 an example of equatorial super¬ 
rotation corresponding to the steady and stable solu¬ 
tions illustrated in Fig. 3, typically at the intersection 
of the red curve and the solid blue line (Box 2). Here, 
u s is the surface angular velocity at the equator for the 
solid planet, u Q is the zonal wind at the atmospheric jet 
center located at altitude with a pressure p 0 , and u Qc 
is the critical zonal wind below which there is no super¬ 
rotation (Box 2). In the analytic model, the jet center 
is characterized by a layer where the solar radiation is 
strongly absorbed: a cloud layer for Venus, a haze layer 
for Titan, and a dust cloud layer for Earth. For Venus, 
the modeled atmospheric angular velocity at the cloud- 
top level is 65 times greater than its solid body rotation. 


The modeled Venus u Q (118 m s _1 ) is very close to the 
measured zonal wind of =115 m s _1 at the 65-km cloud- 
top level. 16 

For Titan’s stratosphere, limited measurements 1 ^ and 
numerical simulations by a 3-D model 13 suggest super¬ 
rotational winds of =100 m s _1 centered at the 1.0-mbar 
(100-Pa) pressure level. My analytic model shows that 
no solution is possible for equatorial superrotation at 
this level, corresponding to the case of nonintersection 
between the blue dashed line and the red curve in Fig. 3. 
Titan is far from the Sun and only intercepts a small por¬ 
tion of the Sun’s energy. However, if the main solar radia¬ 
tion absorption layer in Titan’s stratosphere is lifted from 
1.0 mbar (=185 km) to 0.1 mbar (=288 km), which 
increases the solar heating rate per unit volume by a 
factor 10 because of decreasing air density, an equatorial 
superrotation of 106 m s^ 1 could be maintained at the 
0.1-mbar (10-Pa) level (Table 2). It is hoped that ongoing 
observation by the Cassini-Huygens mission will provide 
suitable data to validate this model prediction. 

The last three rows in Table 2 show the relative con¬ 
tributions to the momentum budget from three terms in 
Eq. 2. Negative values indicate momentum sinks. The 
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BOX 2: STABILITY ANALYSIS OF A BALANCED SUPERROTATIONAL JET AT THE EQUATOR 


The source of zonal momentum is wave pumping, which is 
the dynamical effect of the nonlinear wave/mean flow interac¬ 
tion. It is a nonlinear function of the zonal wind as shown by 
the red curve in Fig. 2. However, the momentum sink can be 
expressed as a linear function of the strength of the jet. The 
strength of the momentum sink relative to the momentum 
source is measured by the slope of the straight lines in Fig. 2. 
The jet center acceleration can be described by the difference 
between the wave pumping and the relative momentum sink, 




— ar ], 


( 1 ) 


where tj is the dimensionless wind at the jet center and the 
two positive coefficients s and a denote the relative strengths 
of the acceleration and momentum sink, respectively. When 
the wave pumping is exactly balanced by the momentum sink, 
zonal wind acceleration vanishes: 


substituting 17 = tj 0 + 17' into Eq. 1 , using Eq. 2 , and keeping 
only the linear terms on 77' since it is a small perturbation, 1 
obtain the following linearized equation governing the evolu¬ 
tion of tj 


£~z~ = ~ <*(4 — ■ (3) 

dt 

When tj 0 = TJq =1/2, the right-hand side of Eq. 3 vanishes. 
Corresponding to the tangent point of the blue and red lines 
in Fig. 3, tj 0c is the critical zonal wind below which there is 
no solution for a steady and stable equatorial superrotation. 
The dimensional values of critical zonal velocity u„ for three 
planets are listed in Table 2. 

The solution of Eq. 3 can be written as 

V = 1 7 t ' = o ex P[ _£_1 “(4 _T ?o 2 ) t ]> W) 


£ ~^ = V ° ie ~ ar, o =0 - 


( 2 ) 


To see whether such a balance leads to a stable or unstable 
solution TJq, I introduce a small perturbation 17 ' that is super¬ 
imposed on TJq to examine how tj' will change with time. By 


where V t = Q denotes the initial perturbation. Equation 4 
indicates that when two solutions exist, say, tj 01 < TJq and 
t}q 2 > tj 0c , from Eq. 2, as shown in Fig. 3, the smaller t)qj corre¬ 
sponds to an unstable solution because an initial perturbation 
V t = Q will grow exponentially, whereas the larger tj 0 , leads to 
a stable solution because of its exponential suppression of an 
initial perturbation. 


Table 2. Equatorial superrotation and momentum budgets for the atmospheres 

Bodies 

of three bodies. 


Venus 

Titan 

Earth 

Rotation rate, (1 (s -1 ) 

3.001 X 10 -7 

4.573 X 10 -6 

7.292 X 10 -5 

Radius, a (m) 

6.052 X 10 6 

2.575 X 10 6 

6.371 X 10 6 

Surface angular velocity, u s (m s -1 ) 

1.82 

11.8 

465.0 

Zonal wind at jet center, u Q (ms -1 ) 

118 

106 

120 

Critical zonal wind, Sq c (ms -1 ) 

14.5 

86.3 

55.7 

Altitude, Zq (km) 

65 

288 

17 

Pressure, p 0 (mbar) 

100 

0.1 

100 

Momentum source, III 

100% 

100% 

100% 

Momentum sink, VI 

-31.5% 

-81.9% 

-69.5% 

Momentum sink, VIII 

-67.2% 

-17.7% 

-29.8% 


momentum source driving equatorial superrotation is 
wave pumping (III) due to the strong absorption of solar 
radiation in an absorbing layer within the atmosphere. 
Wind shear advection by the horizontal branches of the 
Hadley circulation (VIII) is a major momentum sink for 
Venus. For the atmospheres of Titan and Earth, frictional 
drag (VI) mainly balances the momentum source. 


EQUATORIAL SUPERROTATION IN 
EARTH’S LOWER STRATOSPHERE 
AND IMPLICATIONS FOR THE 
EXTINCTION OF THE DINOSAURS 

There is no steady equatorial superrotation in the 
Earth’s atmosphere because it is almost transparent to 
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visible solar radiation. The majority of the solar energy 
is either absorbed at Earth’s surface or reflected back to 
space by clouds. Superrotational winds in the Earth’s 
atmosphere (Table 2) could, however, exist if a layer of 
suspending dust were injected into the stratosphere by 
huge meteorite impacts or by massive volcanic erup¬ 
tions, which would lead to the strong absorption of solar 
radiation by the resulting dust cloud. 

Huge meteorite impacts and massive volcanic erup¬ 
tions have emerged as the dominant theories to explain 
paleo-extinction events, specifically the extinction of 
dinosaurs (e.g., Ref. 18) at the Cretaceous/Tertiary 
(K/T) boundary, about 65 million years ago. Alvarez et 
al. 19 proposed that a 10-km-dia. asteroid struck Earth in 
the Yucatan Peninsula in Mexico and injected a layer of 
debris clouds into the lower stratosphere, blocking sun¬ 
light and causing a rapid cooling that wiped out about 
50% of the biota, including nearly all of the dinosaurs. 
If the extinctions near the K/T boundary were gradual 
or stepwise, a volcanic explanation for the extinctions 
would be more consistent (e.g., Ref. 18). However, in 
either theory, an environment favorable to the survival 
of dinosaurs was disrupted locally and quickly spread 
across the entire planet, resulting in global mass extinc¬ 
tion. Numerical simulations by an aerosol model show 
that the environmental consequences of the impact¬ 
generated dust cloud are sensitive to the horizontal 
transport time, since the effect on dinosaur extinction 
requires the dust cloud to remain in the stratosphere for 
a long enough period to have a global impact. 20 

Based on my analytic model, 1 it is found that an opti¬ 
cally thick dust cloud in the Earth’s lower stratosphere 
centered at the 100-mbar level, such as one created by 
an asteroid impact, could produce an equatorial super¬ 
rotation of -120 m s _1 (Table 2). Such a strong equa¬ 
torial superrotational jet would greatly enhance the 
meridional transport of the dust cloud. Two processes 
associated with equatorial superrotation can cause this 
enhanced transport. First, enhanced heating of the dust 
cloud increases the strength of the radiatively driven 
Hadley circulation and reduces the meridional transport 
time. Second, global-scale waves are expected to be gen¬ 
erated from the jet stream, possibly as a result of instabil¬ 
ity, which can also effectively transport the dust cloud 
poleward. The large meridional spreading of organized 
cloud features shown by the dark horizontal Y-shaped 
4-day waves apparent near Venus’ cloud top (Fig. 1) can 
be considered an example of the presence of global-scale 
waves in a strongly superrotational atmosphere. 8 

In the Earth’s current atmosphere (with no strong 
equatorial superrotational winds), large-scale waves 
excited in the equatorial lower atmosphere are usually 
confined within a narrow band about 15° north and 
south of the equator. These equatorially trapped waves 
play a major role in driving equatorial circulation in the 
stratosphere through wave/mean-flow interaction, but 


have little direct effect on mid-latitude atmospheric cir¬ 
culation. 21 Since the zonal mean winds are driven by 
equatorially trapped waves, their effects are also mainly 
confined to lower latitudes. From the viewpoint of tracer 
transport in the Earth’s lower stratosphere, the tropics 
can be considered isolated from the mid-latitudes. 

Measurements and model simulations suggest that 
the transport time from the equator to high latitudes is 
1 to 2 years in the lower stratosphere. 22,28 On the other 
hand, because of the presence of superrotational wind 
and the associated global-scale wave modes near Venus’ 
cloud-top level, it only takes a few days to transport cloud 
features from low- to high-latitude regions on the planet. 
One would expect similar modes of unstable global-scale 
waves to appear on Earth if the meridional variations 
of strong superrotational wind in the lower stratosphere 
were large enough to cancel out the curvature effect 
that traps large-scale equatorial waves. Therefore, more 
efficient meridional transport of a dust cloud would be 
expected to accelerate the globalization of a substantial 
meteoritic dust cloud, rapidly spreading an environment 
unfavorable to the survival of the dinosaurs on a scale of 
a few days rather than over a few years as in the Earth’s 
present atmosphere. 

CONCLUDING REMARKS 

The comparative approach to studying planetary sci¬ 
ence has received increasing attention and has become 
fruitful as more data from different planets have become 
available. 2/1-26 Because of the extensive diversity of 
planetary bodies, in-depth studies of comparative plan¬ 
etary atmospheres are often conducted on a one-to-one 
comparison basis. In this article, atmospheric equato¬ 
rial superrotation is investigated for Venus, Titan, and 
Earth. By a brief and hierarchical overview of various 
modeling methods, I have demonstrated the feasibility 
and usefulness of the comparative approach of studying 
planetary atmospheres using a common set of physical 
laws. The main results are summarized in Table 2, which 
shows that under appropriate conditions all three bodies 
could maintain stable equatorial superrotational winds 
greater than 100 m s . The single momentum source 
that maintains the superrotation is thermal pumping by 
tidal waves, balanced by the two momentum sinks of fric¬ 
tional drag and the advection of the horizontal shear by 
the Hadley circulation. Furthermore, since I have solved 
the problem analytically, the general parameter depen¬ 
dence of equatorial superrotation on various external 
and internal parameters becomes apparent. The detailed 
derivations and explicit expressions for the solution are 
given in a comprehensive paper. 1 The analytic model 
includes the following important parameters: 

• Sun-planet distance 

• Radius of the planet 

• Rotation rate 
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• Inclination of the equatorial plane 

• Gravity 

• Atmospheric scale height 

• Atmospheric buoyancy frequency 

• Rayleigh friction coefficient 

• Albedo 

• Pressure level at which the absorption of the solar 

radiation occurs 

Studies of the atmospheres of other planets also 
benefit our understanding of the Earth’s atmosphere. I 
have already pointed out one example: the understand¬ 
ing of the Antarctic “ozone hole” associated with the 
catalytic destruction of ozone by chlorine and other 
halogen-containing components benefited in part 
from theoretical studies of the photochemistry of HC1 
in the atmosphere of Venus. 2 One important by-prod¬ 
uct of investigating the equatorial superrotation of a 
slowly rotating planetary atmosphere, such as Venus or 
Titan, is the realization that Earth’s atmosphere could 
also have had equatorial superrotation 65 million 
years ago. Since atmospheric equatorial superrotation 
was first observed on Venus, it has always been asso¬ 
ciated with a slow planetary rotation rate. Therefore, 
the possibility of superrotation and its consequences 
on the faster-rotating Earth have never been consid¬ 
ered before in a similar perspective as that of Venus 
or Titan. My analytic solution to the problem, i.e., 
evaluating all of the momentum terms in Eq. 2, sug¬ 
gests that the magnitude of the planetary rotation rate 
is not critical for the existence of equatorial superrota¬ 
tion on a slowly rotating planet. The dominant cause 
of equatorial superrotation is layered solar heating, the 
rate of which is inversely proportional to the air density 
and increases exponentially with altitude. As a result, 
I am also able to examine the possibility of superrota¬ 
tion and its consequences on Earth. This leads to my 
conclusion that equatorial superrotation in the Earth’s 
lower stratosphere following an asteroid impact could 
have accelerated the globalization of an environment 
unfavorable to the survival of the dinosaurs. 

Planetary atmospheric research is a multidisciplined 
field. Among the three major disciplines—dynamics, 
radiation, and photochemistry—I have focused on atmo¬ 
spheric dynamics in this article and discussed radiative 
energy balance very briefly. The APL Space Department 
has a unique science team whose members are actively 
engaged in all three disciplines. Understanding plane¬ 
tary atmospheres is important both for its own scientific 
merits and for gaining insight into the Earth’s biosphere. 
Looking at atmospheric processes on Earth from the 
perspective of the entire solar system will transform our 
view of atmospheric science. The possibility of equato¬ 
rial superrotation in the Earth’s lower stratosphere and 
its implications for the extinction of dinosaurs discussed 
in this article is one such example. 


I recently encountered the following paragraph that 
describes the possible GCM (general circulation model) 
simulation of equatorial superrotation in the Martian 
atmosphere. 

The grid point dynamic model is based on the LMD ter¬ 
restrial model described by Sadourny and Laval [1984]. The 
discretization scheme conserves both potential enstrophy 
for barotropic nondivergent flows [Sadourny, 1975] and total 
angular momentum for constant surface pressure axisymmet- 
ric flow. The latter property was not included in the original 
terrestrial version but was found to be very important for 
simulating the Martian atmosphere in order to avoid spu¬ 
rious prograde zonal winds in equatorial regions [Hourdin, 
1992b]. 27 

Because the Martian atmosphere is known to have 
suspended dust that absorbs solar radiation and produces 
strong tides, the originally simulated prograde equatorial 
zonal winds in the LMD (Laboratoire de Meteorologie 
Dynamique) Martian model may not be entirely due to 
the model’s spurious momentum source. 
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